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ABSTRACT 

Aims. We describe a new technique for measuring accurate galaxy colours from images taken under different seeing conditions. 
Methods. The method involves two ingredients. First we define the Gaussian-aperture-and-PSF flux, which is the Gaussian-weighted flux a 
galaxy would have if it were observed with a Gaussian PSF. This theoretical aperture flux is independent of the PSF or pixel scale that the 
galaxy was observed with. Second we develop a procedure to measure such a 'GaaP' flux from observed, pixellated images. This involves 
modelling source and PSF as a superposition of orthogonal shapelets. A correction scheme is also described, which approximately corrects for 
any residuals to the shapelet expansions. 

Results. A series of tests on simulated images shows that with this method it is possible to reduce systematic errors in the matched-aperture 
fluxes to a percent, which makes it useful for deriving photometric redshifts from large imaging surveys. 
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1. Introduction 

Accurate measurement of colours of astronomical sources is 
important for many applications. In this paper we describe 
a new technique for colour measurement that is specifically 
geared towards poorly-resolved galaxy images in wide-area op- 
tical surveys. It is often the case that images taken through 
different filters have different image quality (different see- 
ing or pixel scale, for example), and it is important to cor- 
rect the photometry for these effects. An important applica- 
tion is the measurement of photometric redshifts of galax- 
ies (e.g., Baum 1957) from multi-band imaging surveys such 
as SPSS (Adelman-McCarthy et al. H gjffi, FIRES (Labbeet 
al. l2003l Forster Schreib er et al . l2006t . FDF (Heidt et al. l2003l . 
CFHTLS (Brodwin et al. 2006 1, or the up-coming KIDS survey 
on the VST (Kuijken et al., in preparation). 

Ideally, the total flux for each source can be measured in 
each band, after which the colours follow simply. But in prac- 
tice we are limited to measuring aperture fluxes, because of 
pixel noise and confusion by neighbouring sources. Finite aper- 
ture fluxes are affected by smearing by the point spread func- 
tion (PSF) and by pixellation. 

To measure accurate colours does not require good mea- 
surements of the total flux of a source: it is sufficient to deter- 
mine the flux of a particular intrinsic part of the source in dif- 
ferent wavelength bands, and compare those. The aperture need 
not have sharp edges: any weight function will do. The techni- 
cal problem to solve is then how to make sure that the same 
intrinsic part of the source is measured in each band, when all 



that is available is a set of convolved images of the source, each 
with a different PSF. 

An often-used solution is to smooth all images to the same 
effective PSF. This involves convolution of the sharpest images 
in a multi-band set with appropriate convolution kernels, result- 
ing in smoothed images that have the same PSF as the worst- 
seeing image of the set. Matched-aperture photometry on these 
images then yields proper colours. A practical disadvantage of 
this method is that all images need to be brought together in 
order to compute the proper convolution kernels — in case of a 
large multi-band survey this is not a simple task. It is also com- 
putationally expensive to convolve all pixels of an image with 
a spatially variable kernel, particularly if most of these pixels 
will in the end not be part of any source. 

This paper describes a different approach which is consid- 
erably more efficient and no less accurate. Rather than con- 
volving images to each other's PSF, we encode the detected 
sources in each image as a sum of orthogonal basis func- 
tions ('shapelets', Refregier 2003 henceforth R03). The same 
is done for the PSF. From this point on we only manipulate 
the sources, removing the need to convolve all pixels. For each 
source we then construct a 'PSF cleaning' operator with which 
to compute the aperture flux that would be measured with a 
standard (Gaussian) PSF. By standardizing on a 'clean' PSF 
we avoid the need to know the 'dirty' PSFs of all bands at once, 
which means that images can be analysed quite independently. 

This paper is organized as follows. Section |2] defines the 
aperture flux we will use. Section summarizes the shapelets 
formalism, and we show there how to use it to compute our 
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aperture fluxes for shapelet objects. An error analysis is also 
given. In section@]we present the results of many simulations, 
and derive several technical improvements which take into ac- 
count the residuals to the shapelet description for the observed 
sources and PSFs. Tests of the noise propagation and a discus- 
sion of the choice of aperture radius are treated in section [5] 
After a short discussion (Section|6} we summarize the paper in 
Section0 

2. An aperture flux measure independent of PSF 

We define the 'Gaussian-aperture-and-PSF' (GaaP) flux F q of 
aperture radius q as the Gaussian-weighted flux the source 
would have if it were observed with a Gaussian PSF of the 
same width as the weight function: 



ee j dre- r2/2 " 2 j dr'S(r')- 



-(r-r')W 

2nq 2 



(1) 



where S is the intrinsic (pre-PSF, pre-pixelization) light dis- 
tribution of the source. We will show below how we can esti- 
mate F q from observed (PSF-convolved, pixelated) images of 
the source. 

Changing the order of integration in eq.^ an d some simple 
algebra shows that F q is equivalent to the integral 



H 



e- r2/4q2 S(r), 



(2) 



i.e., half of the weighted flux of the intrinsic source S with a 
Gaussian weight function of dispersion 2 l ^ 2 q. 

The GaaP flux for Sersic model light profiles is plotted in 
Fig. [0 and compared with the top-hat aperture flux. As ex- 
pected (because the GaaP flux includes some smoothing), the 
GaaP flux has a flatter curve of growth than the straight aperture 
flux, and is less sensitive to differences in galaxy type. 

3. Shapelets 

In the shapelets formalism (R03) images are expressed as a 
truncated sum of orthonormal Gauss-Hermite terms. In this pa- 
per we follow the notation used in Kuijken ( 2006 henceforth 
K06), where an application of shapelets to weak gravitational 
lensing is described. 

The shapelet basis function of scale B and x,y orders a, b is 



B ab (x,y) = Jf b /T 1 e- Kx - x ^ +(y - y ^ y: ^H a (x/j3)H b (y/B), 



(3) 



where x c and y c are the center of the expansion, B" b is the basis 
function of order (a, b), H" is the Hermite polynomial of order 
a and 



k ab = \l^na\b\2 a+b (4) 

is a normalization constant chosen so that (B ab ) integrates to 
one. 

An image of a source, or a point spread function, can be ex- 
pressed as a sum of such basis functions with coefficients s„t,. 
This can be done in two, not quite equivalent ways: by integra- 
tion or by fitting. Integration is more straightforward in princi- 
ple. The B ab are an orthonormal basis, so the coefficient of each 




q/R. 



Fig. 1. Left: The GaaP aperture flux for different Sersic profiles 
of effective radius 1 , and total flux 1 , as a function of aperture 
radius q. From steepest to flattest, Sersic index 0.5 (Gaussian), 
1 (exponential), 2, 3 and 4 (de Vaucouleur) are plotted. Right: 
the top-hat aperture fluxes within radius q for the same profiles. 



term is simply the integral of the image times this basis func- 
tion. Since the image is pixellated, however, performing this in- 
tegral numerically is fundamentally inaccurate since it needs to 
be approximated by a finite sum. Fitting has the advantage that 
it is well-defined: one minimizes^ 2 over all coefficients, which 
yields a description of an image that fits the fluxes in each im- 
age pixel. Formally, the best-fit coefficients therefore describe 
the pixel-convolved image, which means that when image and 
PSF are both fitted in the same way, their deconvolution should 
be independent of the pixel scale. 1 In this paper we determine 
shapelet coefficients by least-squares fitting to the pixel values. 

Once the shapelet coefficients have been determined they 
can be operated on to calculate the effects of convolution, trans- 
lation, magnification, shear, etc, as well as a range of photo- 
metric quantities. For the present application the relevant oper- 
ations are convolution, and weighted-aperture photometry. 

Convolution of a source by a PSF is a linear operator and 
can be represented by the action of a PSF matrix P. The ele- 
ments of P depend on the shapelet coefficients p„b and scale 
j8p S f of the PSF, on the shapelet scale fi m of the source being 
convolved, and on the shapelet scale B oul of the result. For a 
model source whose shapelet coefficients are m a b, the convolu- 
tion with the PSF can be written as 

1 A third possibility is to fit the image as a superposition of pixel- 
convolved shapelet basis functions. This yields a description of the 
pre-pixellation image. A disadvantage of this approach is that the ba- 
sis functions are no longer orthogonal, implying that the coefficients 
depend on the order at which the basis function set is truncated. 
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The coefficients of the PSF matrix are 

p , (ft- r \ - V rP°°tPs* t P™rP<»«0p»fft» n 

1 fl]<!] a 2 b 2 \Pm >PoutJ — /j^aia 3 a 2 b\b$bi ' 



a 3 b 3 



(5) 



(6) 



a 3 ,b 3 



where the elemental convolution matrix C^f 1 expresses the 
convolution of two one-dimensional shapelets of scales ft\ and 
P2 as a new shapelet series with scale ft^,. A recurrence relation 
for C m „i is given in R03. 

If the image of the source S is expressed as a shapelet se- 
ries, with scale radius ft, then the GaaP flux defined in |3 can 
be calculated analytically for each S = B ab : only those terms 
for which the orders a and b are both even contribute, and in 
this case the x and y integrals in eq.|3 separate as F q = F a q F b q 
with 



-x 2 /4 9 - 



H a (x'/ft)e 



^^2 a a\ft 

which, with the help of the result that, for even n, 



e-' 2 H"(tx)dt = V^-^ttO 2 - 1)" 



12 



(n/2)\ 

(Abramowicz and Stegun iTl)72l eq. 22.13.17) reduces to 

3 2\"/2 



1/4/31/2 



F a = 

q 



(a-l)U 2q-p 



VI + ft 2 /2q 2 " all \2q 2 +ft 2 



(7) 



(8) 



(9) 



(where 0!! = 1!! = landn!! = n(n-2)\\). Hence the Gaussian- 
aperture, Gaussian-PSF flux of an object with shapelet coeffi- 
cients Sab at scale q is simply the sum 



w= Z 



S a b 



q q' 



(10) 



a,b even 



To calculate the GaaP flux of an object S observed with 
PSF P, we therefore 

1. express both as shapelet series s a b and p a b, using scales ft 
and ft ps [ respectively; 

2. construct the PSF matrix P using eq. |6] using the same 
value ft as the input and output shapelet scale; 

3. calculate the shapelet expansion for the deconvolved image 
P- 1 • s; 

4. calculate the GaaP flux using eqs.l9landllOI 



The resulting expression is 

z z ^ 



-] 



b\aibi q q 



(ID 



a\,b\ even fl2,^2 



From this expression it is also straightforward to calculate the 
variance of F q : for background-limited imaging the standard 
errors on the s a b are uncorrected and, because of the normal- 
ization of the basis functions, equal to the standard error cr on 
the flux in a pixel. Hence the variance of the GaaP flux is 



a 2 ,b 2 



Var N=Z Z P ^ F q' F '. 



V- 



\a\,bi even 



(12) 
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Fig. 2. The ratio of the measured GaaP flux to the true value, for 
a series of simulated Sersic profile model galaxies convolved 
with Moffat PSF profiles. Models are identified in Tabled The 
horizontal axis is the aperture radius q, divided by the shapelet 
scale used to describe the observed image. See the text for de- 
tails. Shapelet order 8 was used. (Colour versions of this and 
similar figures indicating the parameter dependence of the re- 
sults are included in the on-line materials.) 

Table 1. PSF and Galaxy profiles used in the tests 
PSF type 






Pseudo-Airy (eq.|23j 


1 


Pseudo-Airy with 10% spikes 


2 


Moffat, m = 2 


3 


Moffat, m = 3 


9 


Moffat, m = 9 (nearly Gaussian) 


Galaxy type 







Gaussian 


1 


Sersic, n = 1 (exponential) 


2 


Sersic, n = 2 


3 


Sersic, n = 3 


4 


Sersic, n = 4 (de Vaucouleur) 


5 


Spiral feq.l22l 



4. Tests and improvements 

To test this procedure, we performed the following simulations. 
First we generated a series of model galaxy images with Sersic 
(1968) profiles exp(-kr 1 ^"). These were then convolved with 
Moffat profile (1 + r 2 /a 2 y m PSFs, and pixellated to yield simu- 
lated images. A corresponding pixellated PSF image was gen- 
erated in the same way. All images were built by monte-carlo 
selection, as described in K06. 

Each pixelated image and its PSF was then fitted with a 
shapelet series to order 8. Following K06, we used a scale ra- 
dius ft that is 1.3 times the dispersion of the best-fit Gaussian: 
this choice favours a more accurate description of the outer 
wings of the images, and helps to suppress any sub-pixel scale 
freedom the shapelet fit might otherwise have. At this point it 
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Fig. 3. As fig. |2 but now including the additive correction for 
residual pixel flux to the observed sources' shapelet expan- 
sions. 



Fig. 4. As fig. but now furthermore including the multiplica- 
tive correction for residual pixel flux to the PSFs' shapelet ex- 
pansions. 



is assumed that the correct background level is known; this as- 
sumption is addressed in 5j6] 

The GaaP flux is then calculated using eq. and com- 
pared with the exact GaaP flux, as calculated directly from eq.|2] 
and plotted in Fig.^ We repeated this procedure for a range of 
aperture radii q for each model, and ran a large set of mod- 
els of different Sersic index for the galaxies (0.5 to 4), Moffat 
index for the PSF (2, 3 and 9), galaxy effective radius (from 
completely unresolved to 6 pixels), PSF FWHM (3 to 6 pix- 
els), and location of the galaxy centre with respect to the pixel 
grid (center and corner of a pixel, middle of an edge). 

The results of this first test are plotted in Fig.|2l which show 
that, provided the Gaussian aperture radius q is within a factor 
of two of the shapelet scale of the observed source, the pro- 
cedure yields GaaP fluxes that are accurate to 10% or better. 
Going to shapelet order 12 improves the results somewhat (not 
shown), but not spectacularly. 

While useful, these residuals are still troublingly large. We 
now describe two simple refinements which take account of the 
residuals to the shapelet fits to observed source and PSF. 



(q 2 - g 2 , f )^ 2 , and computing the Gaussian-weighted aperture 



'psf 

flux of this convolution 

-(r-r') 2 /(2 9 2 -2£ sf ) 



fdr f dr'R(r'Y- 

J J 2n(q 2 -g 2 ) 



(14) 



which, after changing the order of integration can be written 
(as a generalization of eq.13 as the aperture integral 



I 



drR(r) q , e ->W-2&) 



Thus we add 



F res - R U 



1 f ,-(^y 2 i )/^i 2 -Kf ) 



(15) 



(16) 



to the result obtained from eq.^2 F° r the simulated Sersic pro- 
files this term can be of order 10% of F q , particularly for very 
peaked images. The effect of incorporating this correction can 
be seen in Fig. [5] 



4.1. Residuals to the observed source expansion 

The first improvement is to take into account the residual 
galaxy flux Rij at pixel (/, j) to the shapelet expansion: 



S ob Xxi, yi )- 



•2* < *2* a Wy). 



(13) 



u.h 



Thus far we have simply ignored this flux, but a better approxi- 
mation is to include its Gaussian aperture flux (albeit with only 
an approximate PSF correction) in the total. Approximating the 
PSF by a Gaussian of dispersion g ps f , we can estimate the resid- 
ual GaaP flux by convolving Rn with a Gaussian of dispersion 



4.2. Residuals to the PSF model 

The second effect is caused by residuals to the PSF's shapelet 
expansion, which can introduce a multiplicative scaling error 
during the deconvolution. (Essentially, if some of the flux of 
the PSF is missed in the shapelet model then deconvolution by 
this underestimated PSF will result in an overestimated galaxy 
flux.) To estimate the effect of this missed PSF flux, we com- 
pute the correction factor for a Gaussian galaxy. 

For a Gaussian galaxy, of dispersion g and unit total flux, 
the GaaP flux F q is equal to (see eq.|2j 



dr-- e -'-/ 4 ?- 

2 2ng 2 2q 2 



<t 



(17) 
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Fig. 5. As Fig.@] but for a two-armed spiral galaxy 



However, if a part P K s(x,y) of the PSF P was not included in 
the shapelet model, then to first order in these residuals F q will 
have been overestimated by an amount 



J dr [P K 



P ® G(q) ® G(g)] < 



-r 2 /2q 2 



(18) 



Here G(cr) is a Gaussian of dispersion <x, ® denotes convolution 
and deconvolution. We can estimate this integral by replac- 
ing the PSF P by a Gaussian whose dispersion g ps f is chosen to 
be a reasonable match to the PSF size. Combining all Gaussian 
(de)convolutions reduces ea.H8lto 



dr 



>G( > /9 2 -4f + « a )' 



(19) 



or, after the by now customary swapping of the order of inte- 
gration, 



<t 



2q 2 +g 2 -g 



f 



drP^e" 2 ^ 2 * 2 - 2 ^. 



(20) 



The ratio of eqs.[3U]and[2]is then an estimate of the fractional 
error that results from the incomplete PSF model. Correcting 
for this error is accomplished by dividing F q +F les (from eas.ll II 
and!16> by 



/psf = 1 + 



2q 2 +, 



2q 2 + , 



?psf 



I 



drP„ 



-r 2 l(Aq 2 +2 g 2 -2g 2 ) 



(21) 



In practice we estimate g, the intrinsic Gaussian size of the 
source, from the best-fit Gaussian radius of the observed source 
g ohs as g 2 * g 2 obs - g 2 psf . 

In figure 0] we show the result of implementing these two 
corrections: the GaaP flux is now systematically correct to 
better than 2%, over a wide range of Sersic index and PSF 
size, and for GaaP aperture radii between half and twice the 
shapelet scale of the observed images. The error is below 1% 
for 0.6 < qlfi < 1.6, with the worst errors arising for PSFs and 
galaxies that are only a few pixels across. 



Fig. 6. As Fig. 0] but for PSFs with diffraction rings. PSFs 
and 1 are the same except that PSF 1 includes four diffraction 
spikes which together contain 10% of the flux. 



l.i 



1.05 



. Galaxy Gl.4,5 PSF P0-P9 (elliptical) 



V; 
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0.95 
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Fig. 7. As Fig. El but for elliptical PSFs of axis ratio 3:2. PSF 
and galaxy types have the same meaning as before. 

4.3. More complex morphologies 

Having established that the method works for circular, regular 
galaxies, we extended the tests to more complicated shapes. In 
Fig. 13 we show the F q measurements for a simulated logarith- 
mic spiral galaxy, for the same PSFs as above; the model used 
is obtained by replacing the random azimuths 9 of stars in a 
circular exponential disk model by 



9 - - cos (2(9 - 



lnr)). 



(22) 



5 is chosen to give maximally sharp spiral 



(The amplitude 
arms.) 

In Fig.[6]we show the results for all galaxy models and with 
two PSFs that mimick diffraction-limited PSFs. PSF type is 
defined as 



PSF0(r) = 



sin r/r 

2 _/_3 



r< 1; 
r > 1. 



(23) 
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and PSF type 1 is the same except that 10% of the flux is put 
in four thin linear 'diffraction spikes' along which the density 
varies as PSF0(r/2^ 2 ). Even here the GaaP flux is recovered to 
better than 1% for aperture q between 0.7 and 1.5 times /3 b s ■ 
In the case where the PSF FWHM is 2 pixels (not shown) the 
differences can rise to 5%. 

Finally, a series of simulations were analysed for circular 
galaxies that were convolved with elliptical PSFs (axis ratio 
3:2). Also here, the results are satisfactory (see Fig.0. 

We therefore conclude that the method works. 

More details of the relation between the residuals and the 
model parameters are given in colour figures included with 
the on-line material. We close this section with plots (Fig. [8] 
and [9ji of the PSF and galaxy models used in the simula- 
tions, and remark that even though the shapelet formalism is 
built around a Gaussian 'parent' flux distribution, it can be 
used succesfully to do photometry on much more extended 
sources. Nevertheless, it may be profitable to consider expan- 
sions around more galaxy-like parent functions in the future. 

5. Noise 

A similar set of tests can be used to investigate the effect of 
noise on the GaaP fluxes: we simply add a constant Gaussian 
noise level to all pixels in simulated galaxy images, as appro- 
priate for background-limited images, and repeat the analysis 
many times. In Fig. ^| we show the comparison between the 
errors prediction from eq.^1 an d the scatter from a large num- 
ber of realizations, using the models of Figs.|2]to[6] 

For adequately resolved (FWHM 3 pixels or greater), or 
more Moffat-like, PSFs, the theoretical error bar adequately 



Fig. 9. Contour plots of the models used in the simulations. Top 
two rows: PSF models 0-9; bottom two rows: Galaxy models 
0-5. Contours levels differ by a factor of two. The box sizes are 
six times the FWHM for the PSFs and twelve times the effec- 
tive radius for the galaxies. See TableE]for detailed descriptions 
of these models. 



predicts the scatter in the measurements, and any bias in the 
results is neglegibly small. 

Deconvolution will always amplify photon noise, so it is 
important to try to choose an optimal aperture radius. This is 
addressed in Fig. \^\ which shows how the S/N of the Gaap 
fluxes and of their ratio varies with aperture radius. In this cal- 
culation we assume that a galaxy is observed with PSFs of two 
different sizes, so that a given choice of aperture radius q re- 
quires a different degree of (de)convolution. The ratio of the 
two GaaP fluxes at a given q represents a colour measurement, 
with the S/N ratios on the fluxes adding in inverse quadrature. 
As the figure shows, there is an optimal aperture to use, but 
the maximum is quite broad, so that useful colours can be mea- 
sured for a wide range of q. It certainly is not optimal to smooth 
the best-seeing image to the worst one. 

6. Discussion 

We have shown that the proposed technique is an effective 
way to generate well-defined aperture fluxes, with accurate er- 
ror estimates, that are independent of PSF and pixel scale of 
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Fig. 10. Top: The scatter in measured GaaP fluxes from sim- 
ulated images, ^jVai{F q ) I F q , versus the theoretical prediction 
of eq.^] Each point represents a set of 100 noisy realizations 
of a particular combination of PSF and Galaxy models from 
Tabled with aperture radius q between 0.6 and 1.4 times /? b s . 
The scatter is adequately predicted for signal-to-noise ratios 
above about three. Bottom: the average GaaP flux for each set 
of 100 realizations, versus the rms scatter, showing that there is 
no noise-related bias in the results. 



the observations. In all cases with PSFs that are well-resolved 
(FWHM at least 3 pixels) systematic errors are below a percent. 

The advantage of this approach is that it allows multi- 
colour photometry to be done on a distributed dataset. All that 
is needed is a set of object coordinates and aperture radii, and 



Fig. 11. The effect of aperture radius on S/N of the GaaP fluxes 
and colours. The thin solid curves show the S/N of the GaaP 
flux of an exponential disk, R e — 2 pixels, convolved with 
Moffat j3 m = 3 PSFs of FWHM 3 and 6 pixels (top and mid- 
dle curves — tick marks indicate the shapelet scale radius used 
to model the sources). The heavy line shows the S/N of the 
ratio of the fluxes. For this calculation is was assumed that 
both images have the same pixel noise of 1% of the total flux 
of the galaxy. The dashed line shows the same result for a de 
Vaucouleur model galaxy. 



PSF models for each image: then tables of GaaP fluxes can be 
extracted from each separate image and simply merged into a 
multi-colour catalogue. The ratios of the GaaP fluxes of aper- 
ture radius q are then true integrated colours, corresponding 
to the ratio of the Gaussian-weighted fluxes of the intrinsic 
sources with the same weight function exp(-r 2 /4g 2 ). It should 
be stressed that the GaaP fluxes are not total fluxes. 

This method does not absolve one from knowing the cor- 
rect background level and zeropoint of each image. This can 
be difficult with extended galaxies and PSFs, and any flux that 
is missed in either will have repercussions for the measured 
aperture fluxes. Provided the images are well flat-fielded, it is 
possible to obtain decent local background estimates in annuli 
around each source or star. This will include some source flux 
from the wings of the flux profile, but for sufficiently large an- 
nuli the effect is small. We have tested the effect with the simu- 
lations described in ^4] by subtracting a background level (de- 
termined as the median flux in the pixels between 8/3 and 1 1/3 
from the centre of the source) before performing the shapelet 
decomposition. The resulting GaaP fluxes, including all correc- 
tions for pixel residuals, are shown in Fig.^] Not surprisingly, 
larger apertures are more affected by background level errors, 
particularly in the PSF. Overall, while the effect is visible, it 
is not catastrophic: moreover to first approximation the same 
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Fig. 12. The same simulation as in Fig.0] except that before the 
shapelet fitting a locally-fitted background level, determined 
from an annulus between radius 8yS and 11/3 is subtracted. 
Particularly the 'wingy' PSFs 0, 1 and 2 are affected by this 
background error. 



amount of 'wing' will be subtracted in different bands so that 
the overall colours should be affected even less. 

We have argued that our procedure is advantageous in the 
context of large multi-band imaging surveys. The alternative 
approach of convolving entire images to a common PSF, fol- 
lowed by matched-aperture photometry, is less practical in that 
it requires much more information exchange between the dif- 
ferent images (it may also be more computer-intensive as it re- 
quires all pixels to be convolved). In the context of huge multi- 
band surveys, in which different bands may be observed with 
different instruments or processed in different places, the logis- 
tics of these exchanges are complex. 

A different, catalogue-based solution to the problem of de- 
riving accurate colours would be to generate model total fluxes 
for all sources, e.g. by fitting detected sources to combinations 
of PSF-convolved Sersic profiles, as done for the Sloan Digital 
Sky Survey (e.g., Adelman-McCarthy et al. 2006 1. This can 
also be done image by image, and also yields fluxes that are 
PSF-independent, but the results will only be valid if the model 
is appropriate. Few-parameter models are a priori unlikely to 
yield percent-level accurate total fluxes, and any colour gradi- 
ents in the source will give rise to wavelength-dependent sys- 
tematic errors to such a simple fit. By contrast, the method 
proposed here makes virtually no assumptions on the galaxy 
shapes. 

The technique described in this paper can be used to com- 
pare any set of images of a piece of sky: while we have concen- 
trated on colour measurements, another application would be to 
compare images taken at different epochs to measure time vari- 
ability. Most time-varying sources are point-like, in which case 
the more straightforward approach of comparing PSF-fitting 
fluxes might be more appropriate, but a possible application 
would be variability of embedded point sources such as galac- 
tic nuclei or distant supernovae. 



7. Summary 

We have defined a PSF-independent aperture flux, the GaaP 
(Gaussian-aperture and PSF) flux, and showed how it can be 
extracted from astronomical images provided the point spread 
function is known. An extensive series of tests on simulated im- 
ages demonstrates that the level of systematic residual is at or 
below the 1 % level, for aperture sizes that are up to a factor of 
two larger or smaller than the observed source. Implementation 
of these ideas into robust data processing pipelines will allow 
accurate colours to be extracted from large wide-field imaging 
surveys, virtually free from residual PSF and pixel scale effects. 
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Appendix A: Colour-coded residual plots 

This Appendix shows the relation between the residuals in F q 
measurements and the model parameters. The data plotted are 
identical to those of Figs. 13 to but the points have been 
colour-coded according to different model parameters. This 
helps to identify the most discrepent cases. 

As a key to the figures, a summary of all PSF and galaxy 
models tested is given in Tabled 
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Fig.A.l. The results of the simulations plotted in Fig.|2] shown here as four panels. Each panel shows the same results, but 
colour-coded according to a different input parameter: top left, Galaxy type; top right, Galaxy effective radius; bottom left, PSF 
type; bottom right, PSF FWHM. The colour used to plot the results for each parameter value is shown in each panel, and the 
model types are identified in Tabled 
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Fig. A.2. As Fig. IA.ll but for the results of the simulations plotted in Fig. [3] 
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Fig. A.3. As Fig. IA.ll but for the results of the simulations plotted in Fig.|4] 
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Fig. A.4. As Fig. IA.ll but for the results of the simulations plotted in Fig.[5](spiral galaxy model). 




Fig. A.5. As Fig. IA.ll but for the results of the simulations plotted in Fig.[6](Airy function-like PSFs). 
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Fig. A.6. As Fig. IA.ll but for the results of the simulations plotted in Fig. (elliptical PSFs). 




Fig. A.7. As Fig. IA.ll but for the results of the simulations plotted in Fig. [21 (locally fitted background subtracted from source 
andPSF). 



